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The effect of impurities situated at different distances from a two-dimensional electron gas on the density 
of states in a strong magnetic field is analyzed. Based on the exact result of Brezin, Gross, and Itzykson, we 
calculate the density of states in the whole energy range, assuming the Poisson distribution of impurities in 
the bulk. It is shown that in the case of small impurity concentration the density of states is qualitatively 
different from the model case when all impurities are located in the plane of the two-dimensional electron gas. 
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1. Introduction. Two-dimensional electrons in a 
quantizing magnetic field H has been attracting much 
attention £Q especially since the discovery of the quan- 
tum Hall effect The properties of two-dimensional 
electrons in the magnetic field are affected by the pres- 
ence of electron-electron interaction as well as by impu- 
rities. Investigation of the density of states as a function 
of the magnetic field and filling fraction allows us to esti- 
mate the inhomogeneities caused by impurities in exper- 
imental samples [3]. Although the electron-electron in- 
teraction should usually be taken into account, the ques- 
tion about the density of states in the simplest model 
of noninteracting electrons is also rather interesting. 

In the absence of interaction, impurities near a 
two-dimensional electron gas (2DEG) provide the only 
mechanism for broadening of Landau levels. In a weak 
magnetic field the large number of Landau levels, N 3> 
1, are filled. One can therefore use the self-consistent 
Born approximation that is justified by the small pa- 
rameter \nN/N <C 1. It results in the well-known semi- 
circle shape for the density of states 4 . Beyond the 
self-consistent Born approximation one can find the ex- 
ponentially small tails in the density of states [5] . 

In the opposite limit of a strong magnetic field only 
the lowest Landau level is partially occupied. In this 
case one can neglect the influence of the other empty 
Landau levels assuming u>h S> T, t -1 . Here luh = 
eH/m denotes the cyclotron frequency with e and m be- 
ing the electron charge and mass respectively, T stands 
for the temperature, and r is the elastic collision time. 
The density of states on the lowest Landau level strongly 
depends on the statistic properties of the random po- 
tential created by impurities and on the value of the 
dimensionless parameter ns/n^, where Ul = l/(27r/jy) 
with the magnetic field length Ih = 1/ and ns 
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stands for the two-dimensional impurity density. For 
the white- noise distribution of the random potential the 
density of states was found exactly by Wegner |Q ■ For 
arbitrary statistics of the random potential the density 
of states was obtained exactly in a beautiful paper by 
Brezin, Gross, and Itzykson 7 . If the number of impu- 
rities is less than the number of states on the Landau 
level, ns < n^, the Landau level remains partially de- 
generate. In the opposite case ns > n^ the presence 
of impurities leads to complete lifting the degeneracy of 
the Landau level 017]. 

In experimental samples the impurities can be found 
rather far from the 2DEG [HE]. In such a situation 
the two-dimensional electron system is subject to the 
three-dimensional random potential. It means that an 
electron localized at the hetero junction feels impurities 
situated at distances much larger than the width zq 
of the 2DEG. This situation was considered recently 
by Dyugaev, Grigor'ev, and Ovchinnikov [Hj. Within 
the lowest order of the perturbation theory in the con- 
centration rtj mp of three-dimensional scatterers, they 
have calculated the density of states D(E) in the limit 
when the multiple scattering on the same impurity pro- 
vides the main contribution. Assuming exponential de- 
cay of the wave function in the transverse direction, 
<p 2 (z) cx exp(— z/zq), they obtained a universal regime 
where D(E) = n- lmp Zo/E, and energy E is measured 
from the unperturbed Landau level. Being bounded 
both from the sides of small and large energies by many- 
impurity effects, this interval contains most of the states 
of the unperturbed Landau level. Though the analysis 
of Ref. 9 holds for an arbitrary Landau level, it cannot 
be generalized to the limits of small and large energies 
where a nonperturbative treatment of impurity scatter- 
ing is required. 
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The main objective of the present letter is to present 
the full analysis of the effect of far impurities on the 
density of states of a two-dimensional electron gas in 
a strong magnetic field. Employing the remarkable re- 
sult of Brezin, Gross, and Itzykson [7] we calculate the 
broadening of the lowest Landau level by the three- 
dimensional short-range impurities with the Poisson dis- 
tribution in the bulk. 

2. Results. Usually impurities occupy rather large 
volume near a two-dimensional electron gas and, conse- 
quently, their number exceeds the number of states on 
the Landau level, N- lmp n^S, with S being the area 
of two-dimensional electron system. Therefore, the de- 
generacy of the Landau level is removed completely by 
impurities (7JE1- The behavior of the density of states 
is determined by the new dimensionless parameter 
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that will be referred to as impurity concentration. Here 
the three-dimensional impurity density and z 
stands for the spatial extent of the electron wave func- 
tion in the direction perpendicular to the 2DEG, explic- 
itly defined in Eq. (gj. 

In experiments there is usually a small amount of im- 
purities in a layer of width zq near the two-dimensional 
electron gas i-e. impurity concentration is small, 

/< 1. In this case we obtain the following density of 
states at the lowest Landau level as a function of the 
deviation E from the unperturbed level ionj2: 



D(E) 



where 
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On deriving Eq. @ we have assumed that the wave 
function ip(z) decays in the transverse direction as 
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ip {z) exp 
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Z > ZQ, 
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with zq defining the width of the 2DEG, and A being a 
constant of order 1. The form Q corresponds to a rect- 
angular well confining potential p^. The energy scale 
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e^A 
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is introduced by impurities, with uq > being the 
strength of the repulsive disorder potential [cf. Eq. (|14f> 
below], and 7 « 0.577 denoting the Euler constant. The 
result for E 3> Eq is governed by the parameters 
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where z\ ~ zq is the width of the wave function as de- 
termined via its fourth moment: 



1 



(p (z)dz. 
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The fact that the density of states vanishes for E < 
is expected since the random potential is purely repul- 
sive. Since lim^oo 5/(£) = 0, the density of states 
also vanishes at the position of the unperturbed Lan- 
dau level, D(0) = 0. 

In the interval < E <C Ene" 1 ^ the density of 
states exhibits the maximum 
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at the exponentially small energy 
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In the region Ene~ x ^ <C E <C E the density of 
states D(E) = fn^/E is linear in impurity concentra- 
tion coinciding with the perturbative result obtained by 
Dyugaev, Grigor'ev and Ovchinnikov It indicates 
that the multiple scattering on the same impurity pro- 
vides the main contribution to the density of states for 
energies Ene-Vf < £ < £ - This energy interval con- 
tains the major part of the states formed from the lowest 
Landau level. 

The result announced in Eq. J2J in the limit < 
E < Ener 1 '! is applicable for fln(E Q /E) - 1 > yX 
cf . Eq. (|31|l . At the border of applicability Eq. J3J) gives 
Sf[fln(E /E)} ~ fEo/E, and thus D(E) merges with 
the universal result at E ^ Ene~ Xl ' . 

In the region of rather large energies E Eq, the 
tail of the density of states is described by the same ex- 
pression as if all impurities were situated in the plane of 
the 2DEG, with the effective two-dimensional parame- 
ters 



(2D) _ uo_ 
*° ~ Zl 



(2D) _ 
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We mention that the tail of the density of states corre- 
sponds to some optimal fluctuation of the random po- 
tential as it happens for the pure two-dimensional prob- 
lem man]. 
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For large impurity concentration, / 3> 1, the Poisson 
distribution can be replaced by the white-noise distri- 
bution of impurities on the plane with the effective pa- 
rameters IjlUf) . The density of states is given therefore 
by the well-known formula |H1[7] 



D{E) 



n L __ w (E- hEl 



where we introduce the function 



(11) 



W(z) = ~^e z 
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The shift of the maximum of D(E) to positive energies 
is related to repulsive character of impurities' potential. 
Eq. describes the density of states for the Poisson 
distribution only approximately since the exact density 
of states should vanish E < 0. However, deviation of 
Eq. from the exact answer is exponentially small 
(e~^ <C 1) for positive E. 

3. The model. The spin-polarized two- 
dimensional electron gas in the presence of the 
random potential V(r, z) and the strong perpendic- 
ular magnetic field H is described by the following 
one-particle Hamiltonian 

H = (V - ieA) 2 + V(r, z) + U coni (r, z). (13) 

2m 

Here A stands for the vector potential, H = rot A and 
U con f(r, z) denotes the confining potential that creates 
the two-dimensional electron gas. We use the units such 
that h = 1 and c = 1. 

We assume that impurities situated near the two- 
dimensional electron gas arc zero-range repulsive (uq > 
0) scatterers producing the random potential 

JV lmp 

F(r) = u <5 (2) (r - T 5 )S{z - z ). (14) 

i=i 

Assuming that the confining potential U con { depends 
only on the z coordinate, we can represent the electron 
wave function as follows 



(15) 



where tp{z) is the ground-state wave function for the 
electron motion in the direction perpendicular to the 
2DEG in the absence of disorder, and ip(r) describes the 
electron motion in the plane of 2DEG. The decomposi- 
tion (|15fl is equivalent to the projection onto the lowest 
level of dimensional quantization and is analogous to 
the projection onto the lowest Landau level states tp(r). 



Since in experiment the energy separation between the 
lowest and the first excited level of dimensional quanti- 
zation is usually larger than the cyclotron gap, the accu- 
racy of projection onto f(z) is higher than the accuracy 
of projection onto the lowest Landau level. With the 
help of the ansatz l|15fl the original three-dimensional 
problem (|13fl reduces to the two-dimensional one with 
the effective two-dimensional random potential 

AW 

V cS (r)=u J2<P 2 ^) si2) (r-r J ). (16) 
i=i 

Thus, the distribution of impurities along the z di- 
rection leads to an additional random distribution of 
the potential strengths uo<p 2 (zj) effectively felt by two- 
dimensional electrons. 

By using the general result of Brezin, Gross, and 
Itzykson [7] for the random potential (|16fl . we obtain 
for the density of states at the lowest Landau level: 



D(£) = ^ImAlnF(£), 

7T OE 



where 



CO / t 

\ 



(17) 



(18) 



The properties of the random potential are encoded in 
the function g(P) which is defined as 



exp^ n L 



0]} 



exp 



d 2 r/3(r)K ff (r) 



(19) 

where the average (• • • ) is with respect to the distribu- 
tion of the random potential V c s (r). 

We assume that the three-dimensional scatterers 
(|14fl with equal strengths uq obey the Poisson statistics, 
being uniformly distributed along the z direction. Then 
averaging over V e s(r) in Eq. (|16|l reduces to integration 
over z coordinate: 



"'imp 



(20) 



On writing Eq. H2U|) we employed the fact that the wave 
function ip(z) vanishes for z < 0. 

4. Evaluation of the density of states. The 
density of states is generally given by the integral repre- 
sentation (JTZJl, CHJ an d P"V However, Eq. cannot 
be calculated in a closed form valid for arbitrary val- 
ues of impurity concentration and energies. Below we 
analyze the most interesting asymptotic cases. 
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First of all, we note that D(E) vanishes for ener- 
gies E < regardless of the form of <f(z). This follows 
from the fact that for E < the function F{E) is purely 
imaginary that can be obtained by performing the Wick 
rotation t — > —it of the integration contour in Eq. (|18fl . 

The density of states can also be easily calculated in 
the limit of either large impurity concentration (/ 3> 1) 
and arbitrary energies, or small impurity concentration 
{f 1) but large energies E ^> E - In both cases 
the integral (|18fl is determined by small values of t that 
allows to expand the function g((3) given by Eq. i|20|) : 



d/3 



''imp 



u t 



2+2 



''imp 



ttgt 
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(21) 



where z\ is defined in Eq. Q. The quadratic term 
in Eq. (|21(l describes Gaussian (white-noise) distribu- 
tion of impurities , whereas the linear term accounts 
for the energy shift due to the nonzero average poten- 
tial of impurities. Employing the result of Ref. 0, we 
arrive at Eq. (|llfl . Using the asymptotic expression 
W{x) ~ 2y/irx 2 e~ x valid at x 1, we obtain the result 
© for / < 1 in the regime E » E . 

The most interesting is the behavior of D(E) in the 
limit of small impurity concentrations, / C 1, and suf- 
ficiently small energies, E <C Eo- In this limit assumed 
hereafter the function F(E) given by Eq. H18|l is deter- 
mined by large values of t that allows to use the asymp- 
totic formula (@J for calculation of <?(/?) in Eq. (|20|l . 
Introducing the dimensionless energy e = E/Eq where 
the energy scale Eo is defined in Eq. (|5j) and rescaling t 
accordingly we rewrite the expression for the density of 
states as 

D(e) = ^lm^-\nF(e), (22) 



where 



F(e) = J dte le '" £t e- fh ^, 







(23) 



(24) 



The function h(t) is positive at the negative part 
of the imaginary axis, t — —ir, having the following 
asymptotic behavior at r ^> 1: 

h(-ir) = X - lnVr) + c Q + h(r), (25) 

where Co is a constant of the order 1, and h{r) decays 
exponentially at large r: 



h(r) = 



dx 



In a; 



1 



(26) 



The In 2 t asymptotics of h(t) is specific to the problem 
with distributed strengths uoip 2 (zj) of impurities and 
asymptotic behavior Q of the wave function ip{z) far 
from the 2DEG, and should be contrasted with the lni 
dependence for the case of the Poisson distribution with 
constant impurities' strengths. For another decay law of 
the wave function, tp 2 (z) ~ exp[— (z/zo) a ], the leading 
asymptotics would be h(t) ~ ln 1+1 ^ Q t. 

The function F(e) in Eq. I)23|) is given by an oscil- 
lating integral. Therefore it is desirable to deform the 
integration contour to get rid of oscillations. However, 
for e > such a deformation in Eq. I|23|) is impossi- 
ble: the first factor prohibits deformation into the lower 
half-plane, whereas the second factor leads to a diver- 
gent integral if deformed into the upper half-plane. This 
complication can be overcome by splitting the integrand 
into two parts singling out the leading log-square asymp- 
totics: 

oo 

F(e) = Ute ielEt -^'^^^[l+{e-^-l)], (27) 



where we have omitted an irrelevant factor e ^ c °. In 
the limit e < 1, the integral with the second term 
( e -/M«) _ l) in the square brackets allows deformation 
of the contour to the negative part of the imaginary axis, 
where the integrand is purely real. It can be shown that 
the resulting contribution can be neglected compared to 
the integral with the first term in the square brackets. 
The latter can be calculated by deforming the integra- 
tion contour to the upper part of the imaginary axis. 
After a proper rescaling of variables one finds: 



R:) / dr<-x]»(-T-^[lnl 



(28) 



where we have again omitted an irrelevant factor ie~ 7 /e. 

Equations (|2*2*|l and l(2*51) give the integral represen- 
tation for the density of states at E Eq. Its behavior 
depends on the value of the parameter 



1 



(29) 



For small £ < 1, i.e., not too close to the unper- 
turbed Landau level [er x ^ <C e 1), one can calculate 
F(e) perturbatively. Expanding Eq. l(2"%|l in / one can 
easily recover the perturbative result of Ref. |H] as well 
as the leading correction to it: 



D(e) 



EnE 



l-2C(3)/ 2 ln 



1 



(30) 



Retaining only the leading term we obtain the result J5J 
in the regime EqC 1 ^ <C E <C Eq. 
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For large £ > 1 corresponding to energies close to 
the unperturbed Landau level, evaluation of Eq. 128(1 is 
subtler. In this case the ratio Im F(e)/ Re F(e) is expo- 
nentially small and special care must be taken in order 
to extract ImF(e). On the other hand, Rei 7 '(e) can 
easily be calculated for £ — 1 ^> y 7 / '. Making the substi- 
tution r = ee p and calculating the resulting Gaussian 
integral over p one finds 



ReF(e) 



T ex P 



1 

V 



(3i) 



To extract lmF(e), we find it convenient to pass to 
another representation for the function F(e). To this 
end we decouple the square term in the exponential of 
Eq. ((28(1 by the Hubbard-Stratonovich transformation, 
and integrating over r we obtain 



F(e) 



r(l + iz) exp 



z ■ , 1 
— 7TZ + iz In — 

2/ s. 



(32) 

This representation in terms of the T-function is suit- 
able for numerical simulation due to rather fast conver- 
gence of the integral, contrary to the initial representa- 
tion (23J. 

To proceed we shift the integration contour to the 
upper part of the complex plane: z — i£ + x, with x be- 
ing the new real integration variable. As soon as £ > 1, 
in doing the contour transformation, we have to cross 
the poles of the T-function at z = ik with integer k > 0. 
As a result we obtain 



F(e) 



1 



Kl 



fc=i 



(~l) fc 
(fc-1)! 



■ exp 



2£fc 



2/ 



exp(-g)^,/)}, 



(33) 



where [£] is an integer part of £, and the function $(£, /) 
is defined as 



<*>(£,/) 



dxT(l-£ + ix)e-™eiq> 



x 



(34) 

An advantage of this representation is that the pole 
contribution in Eq. ((33(1 is purely real and hence Im F(e) 
is determined solely by Im<I>(£, /). Employing the iden- 
tity r(l — r})T(r}) — 7r/ sin(7r?7) with r/ = £~ix we obtain 
for the imaginary part of $(£, /): 



oo 

Im $(£,/) = -2vrRc / da; 



exp 



x 



nz - 1*) 



(35) 



Q 




Fig. The density of states D(E) in units of ul/^Eo 
as a function of In E/Eq for different values of impurity 
concentration /. 

In the limit £ 3> 1, the term ix in the argument of the 
T- function can be taken into account as r(£ — ix) ~ 
r(£)e~ ixln £. Thereby we find the following estimate: 



Im $(£,/) 



r(£) 



■ exp 



^ln 2 £ 
2 * 



(36) 



Though Eq. l(3t)|) is formally derived for £ ^> 1, it can 
also be applied at £ > 1 as well, with the error being 
small in virtue of the inequality / -C 1 . 

Now with the help of Eqs. l|5T)l. (|53j> and iJSSJl we 

obtain for £ — 1 ^> \/J: 



F(e) 



= l+i\ 



2/ 



|ln 2 £ 



ReF(e) 

Finally, using Eq. J22J), we find 



r(£) 



D(e) 



E 



Sf(Z), 



(37) 



(38) 



where S/(£) is defined in Eq. Q. Equation 1(38(1 gives 
the result © in the region < E < S e _1// . 

The whole profile of D(E) for e <C 1 can be obtained 
by numerical evaluation of Eqs. 1(22(1 and 1(32(1 . The den- 
sity of states numerically calculated for several values of 
the impurity concentration / is presented in Fig. ^ 

5. Conclusion. 

In conclusion, we evaluated the density of states of 
a two-dimensional electron gas in the presence of the 
strong magnetic field and impurities. The fact that 
impurities are situated at different distances from two- 
dimensional electron gas leads to the dramatic change 
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of the density of states in the case of small impurity 
concentration compared to the case when all impurities 
are situated at the same distance from the 2DEG. 

Using the exact result of Ref. [7] we obtained the 
density of states in the whole energy range for the case 
of the wave function with the asymptotic behavior 
The density of states vanishes at the position of the 
unperturbed Landau level and has a maximum at an 
exponentially small energy ©. The major part of the 
states are localized by single impurities in accordance 
with findings of Ref. [Sj • 

The functional form of the result will be different for 
asymptotic behavior of ip(z) differing from the simple 
exponential decay However, the qualitative struc- 
ture of the density of states is supposed to be preserved. 
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